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Abstract. The interpretation of the optical Appell transformation, as previously elaborated 
^^ I in relation to the free-space paraxial propagation under both a rectangular and a circular 

^N , cylindrical symmetry, is reviewed. Then, the caloric Appell transformation, well known in 

the theory of heat equation, is shown to be amenable for a similar interpretation involving 
the Laplace transform rather than the Fourier transform, when dealing with the ID heat 
equation. Accordingly, when considering the radial heat equation, suitably defined Hankel- 
type transforms come to be involved in the inherent Appell transformation. The analysis is 
aimed at outlining the link between the Appell transformation and the canonical transforms. 
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1 Introduction 

K^ I The Appell transformation [l] is mainly associated with the heat equation (HE) [U EJ [Sj HJ [U El [7] . 

jT ^ ■ Some recent analyses however have shown its relevance in other contexts as well [HJ [9] . 

V^ . In [9] some issues related to the ID HE have been revisited within the context of the free-space 

C^ I paraxial propagation, formally accounted for by the 2D paraxial wave equation (PWE). Thus, 

C^ ' the Appell transformation has been interpreted in the light of the propagation of given source 

functions, which are in a definite relation with the source functions of the original wavefunctions 
(i.e. solutions of the PWE), relatively to both a rectangular and a circular cylindrical symmetry. 
The analysis has been developed by following the Lie-algebra based approach to evolution 
equations, ruled by Hamiltonian operators underlying a harmonic oscillator-like symmetry alge- 
bra, as originally elaborated in a series of seminal papers by Kalnins, Miller and Boyer [lOllllj. 
5^ I In fact, the ID HE and the 2D PWE as well as the (1+1)D Schrodinger equation (SE) 

can be considered as evolution equations, ruled by Hamiltonian-like operators (not necessarily 
Hermitian) which are quadratic in the inherent canonically conjugate variables. 

The ID HE has been the object of deep investigations PEIEIIIIEIEIEIEIE], with the relevant 
Appell transformation being also embedded in the Lie symmetry group |12 | I13 | [T^. Extensive 
studies have been devoted to the PWE and the SE as well. In [15] and [l6], for instance, the 
aforementioned Lie algebra based method has been systematically applied to the 2D PWE, and 
markedly interpreted in terms of solutions obtained by propagating definite "source functions" , 
identified according to that method as eigenstates of specific operators in the inherent symmetry 
algebra \W \ II H [12]. Consequently, the transformations between wavefunctions can be seen to 
trace back to definite relations between the respective source functions, as explicitly shown 
in [17] through an analysis aimed at characterizing transformations between wavefunctions in 
terms of transformations between the relevant source functions. 



*This paper is a contribution to the Special Issue "Symmetry, Separation, Super-integrability and Special 
Functions (S**)". The full collection is available at, http://www.emis.de/journals/SIGMA/S4. html, 
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In particular, the optical Appell transformation represents the mapping between wavefunc- 
tions resulting from Fourier or Hankel pairs of source functions. This has been shown in [S] by 
proving the correspondence through the Appell transformation of wavefunctions generated by 
eigenstates of operators belonging to the PWE symmetry algebra, which are linked by a Fourier 
or Hankel similarity transformation, according to whether a rectangular or circular cylindrical 
geometry is concerned. As a mere consequence, the fractional Appell transformation has been 
introduced by involving the fractional Fourier or Hankel transform. Thus, a family of symmetry 
transformations for the 2D PWE, parameterized by a continuous parameter, has been identified. 

Although a direct check of the result can easily be carried out, the formal steps followed 
through the analysis, as presented in [9] and reviewed below, have been aimed at favoring 
a "visualization" of the process, in view also of the possibility of suggesting a practical scheme for 
the realization of the optical "Appell transformer" as parallel to the optical Fourier transformer. 

The purpose of the present paper is to show that a similar interpretation of the caloric Appell 
transformation can be elaborated, resorting to the bilateral Laplace transform rather than to 
the Fourier transform, as far as the ID HE is concerned. We will frame the analysis within 
the Lie-algebra based formalism, as it - besides its intrinsic formal and conceptual elegance - 
reveals effective in suggesting relations with other fields or theories. By retracing the same 
formal steps followed when dealing with the PWE, we will prove in fact that the caloric Ap- 
pell transformation manifests the correspondence between temperature functions generated by 
eigenstates of operators belonging to the HE symmetry algebra, which are linked by a Laplace- 
similarity transformation. Then, resorting to the definition of the fractional Laplace transform, 
a fractional caloric Appell transformation will as well be introduced. 

Although, as earlier noted, the caloric Appell transformation has already been embedded in 
the Lie symmetry group [El HSl [T3] , to the author's knowledge, the Laplace-based relation of 
concern in the present context has still not been explicitly discussed in the literature. In our 
opinion, indeed, the parallelism between the HE and the PWE, which stems from the underlying 
algebra, is still not thoroughly exploited, whilst it may yield interesting results. For instance, 
the property, proven in [6], that the Appell transformation for an in general n-dimensional HE is 
essentially the only symmetry transformation (in the sense that every symmetry transformation 
can be obtained by composing Appell transformations with suitable scalings and shifts of all 
the variables) has been extended in [18] to the 2D PWE as well, by resorting to the symplectic 
ray-matrix formalism, usually adopted in paraxial optics. The stated property stems from the 
fact that the generators of the PWE symmetry algebra are Fourier or Hankel-similarity related. 

Likewise, in regard to the radial HE we will prove that the relevant Appell transformation 
involves a similarity transformation by properly defined Hankel-type transforms. 

As a result, one can recognize a direct link between the Appell transformation and the linear 
canonical transforms, since the (fractional) Fourier transform is an important representative of 
the real linear canonical transforms as the (fractional) bilateral Laplace transform is an equally 
important representative of the complex linear canonical transforms. Likewise, as we will see, 
the transforms, which relate to the Appell transformation for the radial HE, can be framed 
within the context of the radial canonical transforms. 

In Section [2] we will firstly review the basics of the symmetry algebra-based method and later 
the results presented in [9], where the optical Appell transformation has been introduced and 
related to the Fourier or Hankel transform according to whether the PWE in rectangular or radial 
coordinate is concerned. In Section[3]the basics of the canonical transforms theory, specifically in 
relation with the bilateral Laplace and Hankel-type transforms, will be reviewed. Then, retracing 
the steps of the analysis presented in Section [21 we will prove in Section 2] that the caloric 
Appell transformation for the ID HE connects temperature functions generated by eigenstates 
of operators in the inherent symmetry-algebra, which are linked through a Laplace-similarity 
transformation. A similar relation will be shown to hold also for the Appell transformation for 



Appell Transformation and Canonical Transforms 



the radial HE, accordingly resorting to suitably defined Hankel-type transforms, explicable as 
"radial-Laplace" -type transforms. Concluding comments will be given in Section O 

2 Optical Appell transformation: 
Fourier and Hankel transforms 

There is a well-known formal analogy between the (2-|-l)D SE, ruling the dynamics in two 
space-dimensions of a quantum particle under the action of some potential, and the 3D PWE, 
describing the paraxial propagation of a monochromatic scalar light-field through some medium, 
under the basic correspondences i — )• z, /i — )■ A/27r, and some others less direct, relating the 
momentum of the particle and the potential to the direction of propagation of the signal and the 
refractive index of the medium. Due to such an analogy, in fact, the Lie algebra-based method, 
as developed in a series of seminal papers by Kalnins, Miller and Boyer [13 [HI 12] in connection 
with the solutions of the SE, has been applied as well to the free-space 2D PWE, both "linear" 
and "radial" [EJUS], on account in particular of the results in [10], pertaining the (1-|-1)D SE 
respectively for a free particle and for a particle under a potential of the type cjx^ . 

The free-space paraxial propagation is accounted for by the (free-space) 3D PWE, conve- 
niently expressed in the normalized form: 



2 f& 



.d d' d 



^(?,??,C) = 0, (2.1) 



the wavefunction v{S,, t], Q being intended to belong to the complex space ^ of locally C°° 
functions of the real unitless variables ^, r], (". The latter represent the transverse and longitu- 
dinal (Cartesian) coordinates x, y, and z scaled, respectively, to some characteristic transverse 
scale wq, and to the relevant confocal parameter b = kwQ, k being the wavenumber of the light. 
Equation (j2.ip admits separable-variable solutions in both Cartesian and polar coordinates, 
which are suitable to describe wavefunctions displaying respectively a rectangular and a circular 
cylindrical symmetry. In both cases one deals with a 2D equation, involving only one transverse 
coordinate, the longitudinal coordinate C playing the role of an evolution variable. 

2.1 2D "linear" paraxial wave equation 



The separability of the solutions of (12. ip in rectangular coordinates manifests in the factorization 
of the 3D wavefunction v{^, r], C) in terms of two 2D wavefunctions u{^, () and w{r], C) as 

each satisfying the 2D PWE in one transverse coordinate: 



2 1 



d d 



n(e,C) = 0. (2.2) 



It is formally similar to the (1-|-1)D SE for the free-particle Hamiltonian H = — {1/2)8"^ /d^"^. 
The relevant evolution operator (i.e. the paraxial propagator) is 

so that solutions to (j2.2p . for given initial conditions ti(^,0) = uo{^), are obtained as [19] 



n(e,C) = [Z^pwe(C)^o](0 = ^^^ I e*^^no(e'X, (2-3) 



oo 



under the minimal assumption that uo{^) tends to zero sufficiently rapidly as ^ — )• ±oo. 
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The above integral is well-known in paraxial optics as the Huygens-Fresnel (or, Collins) 
diffraction integral (specifically, for free-propagation) |2mi21|. 

Equation (|2.3p implicitly assumes the equivalence of the representations of the propaga- 
tor ^pwe(C) ^s a Fresnel transform and as an exponential operator, involving the free-Hamilto- 
nian. Evidently, the former is meaningful only if uq{^) is integrable and the pertinent integral 
converges, whereas the latter requires the implied series of derivatives -qcttt to exist and converge 
to a finite value as well. We will not dwell here on the legitimacy of such an equivalence, but 
we will simply assume to deal with functions for which it holds |19j . 



2.1.1 Solving the PWE: the Kalnins Miller Boyer method 

The symmetry-algebra based approach, as developed in [lOl |TT1 [12], identifies a systematic 
method to solve (j2.2p . based on the spectral decomposition of operators belonging to the inherent 
symmetry algebra Q. This is the semidirect sum of the algebra sZ(2,M) ~ sp(2,R) ~ su(l,l), 
spanned by the operators {K_|_,K3,K_}, with 

1.2 -- i f ^d l\ ^. 1 92 



[K+,i?-] =2iK3, [K±,i<3]=±iK±, (2.4) 

and the Weyl algebra W, spanned by |X, P, l}, such that 

X:=e, P--^^' 1 = 1' [X,P]=iT. (2.5) 

The algebra Q = W (B sl{2,M) is sometimes referred to in the literature as wsl{2,M); W is an 
ideal in u;s/(2,IR), as conveyed by the mixed commutators, 

[K+,X]=0, [K+,P]=iX, [K_,X]=-iP, [K„,P]=0, 

[K3,X] = -^X, [K3,P] = ^P, [K±,3,I] = 0. 

The operators (|2.4p and (|2.5p are self-adjoint when acting on the Hilbert space £^(R) of 
complex-valued Lebesgue square-integrable functions on the real line M with the scalar product 
{f,g) = J_°° f{09*{0dS,, the star denoting complex conjugation. 

As formulated in [TD], the spectral decomposition {fx}^ of any K £ Q, 

K/a(C) = A/A(e), (2.6) 

can be used to construct solutions of the PWE having the /as as initial conditions. 

This can be done by evaluating the Huygens-Presnel integral (j2.3p for each fx, provided it 
be in the domain of K_, thus obtaining the propagated functions vxi(,, C) = ^pwe(C)/a(0- 

Alternatively, one can search for the solutions of the equations 

K{C)vx{C,C) = ^vx{^,C) (2.7) 

which directly follow from (12. 6p on account of fx{C) = ^pwe(C)~^'^a(^> C)i ^-iid hence address the 
initial conditions vx{C,0) = fxiO- The operator K{(), given by 

K(^) :=e-<K-Ke<K_^ (2.8) 

can be understood as a Heisenberg-like operator corresponding to K. In fact, the above map- 
ping, which relates "fixed- location" operators K G ^ to "evolving- location" operators K{() £ Q, 
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is reminiscent of the transformation relating the Heisenberg and Schrodinger pictures of quan- 
tum mechanics (involving £^(M)-wavefunctions]l|. More precisely, the operator K(C) should be 
regarded as the "back-evolving location" form of K. 
Since 

oo 

e^Be"^ = ^(adA)^B/j!, with (adA^B = [A, [ . . . , [A, [A,B]]]] 

i=o 

being a j-fold commutator, the K(C)s are finite linear combinations of the fixed-location ope- 
rators (|2.4p and ()2.5p with coefficients depending on the propagation variable C? which so en- 
ters (j2.7p as a parameter. 

Evidently, the above sketched procedure, based on the eigenvalue equations (|2.6p or (j2.7p . 
is equivalent to a separation of variables; in fact, in [10] a well-defined connection between 
separation of variables and Lie symmetries for (|2.2p has been established (see also |12j). 

In jl5| . the procedure has markedly been interpreted in terms of wavefunctions obtained 
by propagating definite "source functions", identified as eigenstates of specific operators in Q. 
Accordingly, the transformations between wavefunctions trace back to definite relations between 
the respective source functions [TTj. Further cases besides those considered in [10] have been 
analyzed in [15] . Thus, a certain class of solutions of (j2.2p has been obtained by propagating 
the eigenfunctions of the operator K^q = K3 — (2/^o)K+, ^0 being an arbitrary parameter. Such 
solutions depend on three independent parameters and basically comprise a complex quadratic 
exponential modulated by the Weber-Hermite function D_2iX-i/2 of suitable argument. 

The same result was previously deduced in [23] on the basis of an ansatz giving the general 
solution of the 2D PWE in a suitable separable- variable form. A similar ansatz have led the same 
authors to identify general solutions of (j2.ip pertaining to a circular cylindrical symmetry |24] : 
the elliptical cylindrical symmetry has been considered in [25] . 

On the other hand, as shown in [T7], a definite relation between the wavefunctions arising 
from the eigenfunctions of the operators K^q and K3 can be established, the former following 
from the latter under the symmetry transformation produced by the operator exp[(i/^o)K+(C)]- 

2.1.2 Fourier transform and optical Appell transformation 

As proven in [9], the optical Appell transformation connects wavefunctions, whose source func- 
tions are Fourier related. In order to review this result, we consider for every operator K £ Q 

the dual operator K, linked to the former through a similarity transformation by the direct or 

inverse Fourier transform operator J^, i.e. K = J^KJ-"^^ or K = J^~^KJ^, both belonging to Q. 
The Fourier transform, which is well known to operate as 

+00 



_ 1 r+co 

f{x)=[Tf]{x):=^= e-^^^'f{x')dx' 

\' ZTT J_oo 



admits the exponential operator representation in terms of the su{l, 1) generators K_|_ and K_: 

J- = gW4g-i(7r/2)[K_+K+]^ (^2.9) 

the factor e*'^'^ allowing the optical transform to be matched to the mathematical transform. 



^The Heisenberg picture of quantum mechanics amounts to the observable [22) 
6(i)=Wt(f)6w(i), 
where tt{t) is the unitary evolution operator, which obeys the relations W (t) = U^^{t) and U^^(t) = U{-'t). 
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By optical transform we mean the linear canonical transform resulting from the map of the 
symplectic group Sp(2^ M) (~ SL{2, M)) into the metaplectic group Mp{2, M), which turns optical 
ray-tracing matrices M = ((-, ^) in the former into the Collins integral in the latter as |20^ [2T] 

^ 1 f"^°° ( i ^ 

v{iXo)=U^A s^v{C,Ci) = ^i^=^J expl—{Ae-2CC' + De)\v{^',C^)dC'■ (2.10) 

The optical system described by the ray-matrix M is therefore seen as an operator transporting 
the wavefunction on the input plane located at Q to the wavefunction on the output plane at (o- 
The metaplectic group Mp(2,M) provides a double cover of the symplectic group 5p(2,M), to 
which the Ms belong. In fact, although not explicitly displayed by convention, a double sign zt 
is implied in (j2.10|l i. thus allowing the metaplectic images of the Ms to close into a group. 



It may be useful to clarify that throughout the paper \/i will be intended to identify the 
principal square root of the imaginary unit: yi = e™'^, as it will be in general for the square 
root of any (complex) number: ^/a = Y^o|e''^''^'"^'^, with arg(a) G (— vr, vr]. 

The diffraction integral ()2.3p follows from (j2.10p in correspondence to the free-section matrix 

T(C)=(J ^), (2.11) 

i.e. Z?pwe(C) =^(J C). 

In virtue of the disentanglement relation for the su{l, 1) algebra generators 

gi/3[K_+K+] ^ gitan(/3/2)K_gisin(/3)K+gitan(/3/2)K_ ^ gitan(/3/2)K+gisin(/3)K_ gitan(/3/2)K+ ^ ^2.12) 

holding for |/3| < vr, the exponential operator in (|2.9p is factorizable in the two equivalent forms 

g-^W2)[K_+K+] ^ ^-^K.^-^K^^-^K. ^ ^-iK^ ^-^K. ^-^K^ _ (2.13) 

They reproduce the two possible implementations of the Fourier transform by optical elements, 
i.e. the 2 f -system and the Fourier tube, respectively consisting of a single lens (of focal length /) 
placed midway between two reference planes separated by 2/ and of two identical lenses (of 
focal length /) separated by /. Both setups are described by the ray-matrix F(/) = {_iif n); 
in particular, / = 1 in ()2.9p and so in (j2.13p . Unless otherwise specified, we will address to as 
Fourier transformer the optical system described by the ray-matrix F(l) = F = (-^ q). 

Since every K S ^ is a linear combination of the basis operators (j2.4p and (j2.5p . by (j2.13p 

also J^KJ^~^and J-"~^KJ^ are linear combinations of the same operators, and hence belong to Q. 

Then, as parallel to the eigenvalue problem for the operator K, signified by (j2.6p and (j2.7p . 

we consider the eigenvalue problem for the dual operator K, expressed by the equation 

K5A(e) = A5a(C), 
for the fixed-location operator, and by 

hc)wx{i,C) = Au;a(C,C), K(C) = e-<^-Ke<^-, 

for the relevant evolving- location operator. 

Let us suppose that K = F¥^F~^ . Accordingly, we see that 



A double sign ± is accordingly implied also in the matching factor entering (|2.9 



Appell Transformation and Canonical Transforms 7 

and hence the w\s are solutions of ()2.2p obtained by propagating the Fourier transformed eigen- 
functions of K. Precisely, the wxs are obtained from the vxs through the "local" transformation 

by the symmetry operator 

A:=e-'^^-Te'^^-, (2.14) 

addressed to in [9] as Appell transform operator. In the light of (j2.8p . it can be understood as 
a (back) evolving-location Fourier operator: A = J^iC)- It depends on the evolution variable C; 
however, such a dependence will not be explicitly displayed. 

As noted in [9], the Appell operator ()2.14p individualizes the optical ABCD system, which 
arises from the composition of a Fourier transformer embedded between two free-space sections 
of length — (^ and C- It is therefore described by the symplectic real ray-matrix 

according to the 2x2 matrix representation of each relevant optical component (i.e. operator jf|. 
As a more explicit form of the Collins integral (I2.10p , the propagation of a signal through an 
ABCD system manifests in the transformation of the relevant wavefunction as |171 [27] : 

n(e, Co) = ^(A B^uiC, C^) = -^e*^«'u r|, + f V (2.16) 

It signifies that the whole action of the system is decomposed d la Wei-Norman [28] into a se- 
quence of a free-propagation by B/A, a lensing by focal power — C, and a scaling by 1/A. In 
fact, u{^, (i + B/A) = [^pwE(Ci + B /A)uo\{^) represents the propagated form of the (effective or 
fictitious) source function uo(C) with which the wavefunction u{S,,Ci) can be associated. Note 
that, since we are working with unitless variables, the off-diagonal entries B and C (usually 
having the dimensions of length and 1/length) are intended to be normalized to b and 1/5. 
Relation (|2.16p holds for ^4 7^ 0. If yl = 0, according to (j2.10p the wavefunction transforms as 

Wo BMtCi) = -^e'^^'u(^,C^) , (2.17) 



Kc d) y/iB \B , 

describing the effect of an optical Fourier transformer with focal length f = B combined with 
a modulation by the phase factor exp {ij^^'^)- 

By ()2.16p we can express the associated functions wx in terms of the vxs according to Appell's 
prescription: 

wxi^X) = ^e4vx(-^,-^). (2.18) 

In the case of M^-, A = means C = 0, which amounts to D = and B = 1. Hence, 
following (|2.17p one recovers the primary relation between the source functions: wx{^,0) = 

^A(e,o). 

As noted in [9], relation (j2.18p can also be deduced by acting on the initial function fo(C) by 
the operator e~*''^"-F, and hence (apart from ^/i) through (12.16P with d = and (^^ ^) = (~{ q) . 



^For the practical realization of such a system by basic optical elements, as noted above, one may resort 
to free-sections and lenses for both the free-propagation by ( and the Fourier transformer. According to the 
analysis in [21], the propagation by the "negative" distance — ^ may be realized by free-sections and lenses as 
well, specifically by a sequence of three suitably designed lenses separated by free-sections of proper lengths. 
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If in turn one has K = !F~^K!F, then 

involving just the inverse of the AppeU operator (j2.14p . i.e. 

It individuahzes the ABCD system 

thus yielding the explicit transformation of the wavefunctions f as into the associated w\s as 
wx{(, C) = ^e*^7;A (f -^) • (2-20) 



As a conclusion, we may say that in general the optical Appell transformation 

^^,C) = ^e4^;(^±i,-i), (2.21) 

maps solutions of the linear 2D PWE into solutions [9]. It is a symmetry transformation for 
that equation, which traces back to a (direct or inverse) Fourier relation between the source 
functions of the solutions it connects. It manifests the action of the evolving-location Fourier 
operator J-{C) or its inverse. 

Note that the transformation (j2.2ip amounts to a mathematical Fourier relation between the 
source functions of the involved solutions. Indeed, an optical Fourier relation between the source 
functions would demand for the trivial change of the phase factor as e~"^'^ — )• e~™l'^ . 

As a basic example, we may consider the operators P and X = TVT~^ . The eigenfunctions 
of P, i.e. the plane waves /a(0 = exp(iA^)/v2vr, A G M, yield the wavefunctions 

^^A(e,C) = 4=e'"«e-^^''^/2, (2.22) 



2-K 

which display the familiar frequency-chirping factor exp(— iA^C/2)- 

Applying the Appell transformation ()2.18p to ()2.22p . one just obtains the wavefunctions 

"'A(^,C) = ;7=e*^^, (2.23) 

one would obtain by propagating the eigenfunctions g\(i^ = 5{i — A), A G M, of the operator X. 
Vice versa, the wavefunctions ()2.23p turn into the propagating plane waves (I2.22p by ()2.20p . 
Another example of Appell pair of wavefunctions is provided by the two kinds of Airy beams 



V'KM(e,C, A) = -i.e'^^-^2c-^+2c; Ai (^^ - _ - A 



4 



V'BB(e, C, A) = e-HT2--+2'.; Ai U - ^ - A 



Ai denoting the Airy function of the first kind [29]. As seen in [30], they are respectively obtained 
by propagating the eigenfunctions of the operators V = 2K+ + P and X = 2K_ + X = TVT~^ . 
Accordingly, the former originate from the cubic phase ^km(C) 0, A) = exp(iA^— i^'^/3)/\/27r [TO] . 
as the latter from an Airy pattern: '0bb('?,O,A) = Ai(,^ — A) [31]; in both cases, A signifies the 
eigenvalue of the involved operator "P or Af . 

It is easily verified that by (J2.18P V'km — ^ V'bb, and vice versa by (|2.20p V'bb — ^ V'km- 
In the literature, the V'bbs are addressed to as Airy beams. Along with other similar 
Airy-related solutions of the PWE, they have attracted a great deal of interest (both ana- 
lytical/numerical and experimental), due to their non-spreading nature and their tendency to 
freely accelerate; some pertinent titles arein[3n[32|[33][3^[35l[36][371[38l[39lH0lliT]. 
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2.1.3 Fractional Appell transformation 

Conforming to the relevant generator K„ + K4., the Fourier transform T can be understood 
as the evolution operator describing the dynamics of a quantum attractive oscillator (optically 
conveyed by a quadratic-index focusing medium) observed at "time" r = vr/2, the motion in the 
optical phase-planci resulting in a clockwise rotation by 7r/2. Thereby, according to ()2.14p . in 
the optical phase-plane the Appell transform amounts to a 7r/2-clockwise rotation preceded and 
followed respectively by a negative and positive ^-shear. 

As is well known, the Fourier transform J^ is a specific determination of the more general 
fractional Fourier transform J^", signified by the integral [42 ^ 



{T^^\{x) := JLJ^S^ I "cirE'e^^(""™^'^+"'^°^'^-2""'V(^'), (2-24) 

the order a specifying the angle (j) = aiT/2, or equivalently by the operator 

^a _ gia7r/4g-Ja(7r/2)[K_+K+] _ l^\a 

Again, the factor e*'^'^' ^ relates the mathematical and optical transforms, the latter being yielded 
by the Collins integral ()2.10p in correspondence to the pertinent ray-matrix F" = (_^° ^ a **™s (a) • 
Both expressions reveal the periodicity 

^Q+4j^^a^ j = 0,±l,±2,..., 

so that the range of a can be limited to the interval (—2,2]. Evidently, the ordinary trans- 
form ()2.9p is recovered with q = 1, i.e. = 7r/2: F^ = F. Also, J^ = I. 

Paralleling T, J-"" is interpreted as the evolution operator associated with the harmonic oscil- 
lator dynamics monitored continuously at "time" r = a7r/2 = (p. Accordingly, the corresponding 
motion in the optical phase-plane amounts to a clockwise rotation by (f). 

As conveyed by (j2.12p . J-"" admits two operator factorizations, which reflect its optical realiza- 
tions as "fractional" versions of the implementations of the ordinary transform. One has, in 
fact, ^" = e^<^/2e-^t'^"('?^/2)K-e-*«in('!^)K+e-'tan('?^/2)K-^for the single lens realization (Lohmann 
type I setup [l^) and .F" = e»'?^/2e-*t^'^W2)K+e-»sin(?i)K_g-itan(</,/2)K+ for the two-lens realization 
(Lohmann type II setup [S]). The former comprises a thin lens of focal length 1/ sm.[(j)) placed 
midway between two reference planes spaced by 2tan(0/2), whereas the latter realizes the dual 
configuration, composed by two thin lenses of focal length co\,{(j)/2) separated by sva.[(j)). 

It has been suggested in [9] to consider the fractional Appell transform operator ^", which 
then, as fractional version of the operator A, has been written in the form 

^a ._ g-iCK_^agiCK- = ^"(^)^ (2.25) 

whose dependence on C, will not explicitly displayed. The ordinary operator (j2.14p corresponds 
to a = 1, whilst the inverse is (^")^^ = At"^, since (J^")""^ = F~'^ . 



*Let us recall that the phase-plane (for systems with one degree of freedom) is understood as the Cartesian 
plane formed by the relevant canonically conjugate variables {q,p), obeying the Poisson-bracket relation {q,p} ~ 1. 

The light-ray coordinates {q,p) of geometrical optics are canonically conjugate variables. They span the optical 
phase-plane, where light rays are represented by points, and accordingly the ray propagation through optical sys- 
tems corresponds to the "motion" of the relevant representative point. For instance, free-sections and thin lenses 
produce under paraxial propagation q and p-shears in the phase-plane, i.e. translations of the ray representative 
point respectively in the q and p direction. 

Note that, in accord with the normalization of the space variables here adopted, the ray- variables {q,p) become 
(^) = {l\^"), the ray optical momentum p, i.e. the angle relative to the z-axis, being scaled to the natural far-field 
divergence i?oo = l/kwo, associated with a beam having wo as a characteristic width. 
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The operator (|2.25p individualizes the ABCD system, described by the matrix 
'l C\ ( coscj) sin(/)\ A — C\ /cos (/) — C sin (/) (l+^^)sin( 



-^" \0 ij \— sin (j) cos (j)J \0 1 J \ — sin (j) cos 4> + C sin ( 

which reproduces the matrices ()2.15p for (p = 7r/2 and ()2.19p for (p = —tt/2. 

Then, according to (|2.15p and (|2.17p . we see that the fractional Appell transformation 
amounts to the replacement of the wavefunction v{^, C) by the wavefunction w{S,, C) given by [^ 

fc A^ e^"^/^ -i e^,^in0 f g sm0 + Ccos^\ 

V COS (p — C sm (p \ cos — C sm cos — C sm (p / 

as far as cos 4> — C sin 7^ 0, whilst in the case when cos (j) — C sin </> = (i.e. C = cot 0, which can 
occur for G (— vr, — 7r/2] U [0,7r/2]), it acts as 

u;(^,C) = ^====e^«^ -— I— -,C for C = cot </.. 

The fractional Appell transformation identifies a family of transformations parameterized by 
the continuous parameter (j) S (— 7r,7r]. It would allow one to pass, for instance, from the V'KMS 
to the ^bbs through a continuum of "states" filling up the direct transformation V'KM — ^ V'BB, 
considered before, occurring in fact for (j) = 'k 12 [9]. 

We conclude by recalling that, as a consequence of the properties of the fractional Fourier 
transform, the fractional Appell operator satisfies the basic properties: 

(i) J^ is continuous for all values of the order parameter a, 

(ii) J^ obeys the group property, so that composing two operators of order a^ and a^ yields 
the operator of order a^ +02- 

t/-^ i/~v ■J'\. — ■J~\. »Al , 

(zii) J^ reduces to the ordinary operator for a = 1: A^ = A, and the identity operator for 
a = 0: A^ =1 

The property (i) follows from the continuity of the Collins integral (j2.10p on account of that 
of the ray-matrix entries as functions of (p (or, a). Property (ii) ensues from the fact that 
the composition of the propagation integrals ()2.10p reflects that of the ray-matrices, and hence 
^«i^a2 ^ M_^,jM^,2 = M_j„j+,2 o ^"i+"2 in accord with the addition formulae of the 
circular functions. Property (iii) is implied by that J-^ = J- and J-^ = I. In addition, as for J^", 
the group property (ii) allows A" to be mathematically understood as the a-th power of A. 

Evidently, the expressions (I2.2ip and ()2.26p conform to the well-known rule, which, according 
to standard results from Lie theory, conveys the action on functions in J of any operator in the 
PWE symmetry group spanned by |K_|_,K3,K_} |10tll2j. As earlier noted, the analysis above 
favors the "visualization" of the Appell transformation in optical terms, and makes also explicit 
the connection between some given transformations, which so appear as its "fractional powers". 

2.2 2D "radial" paraxial wave equation 

As we know, the 3D PWE in (normalized) circular cylindrical coordinates (p, 99, ("), 



d d' Id Id 



2?Tr: + Tr^ + -7^ + 



dC dp^ p dp p^ dip^ 



u{p,^X) = 0, 
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with p = vs+^ and ip = arctan(ry/^), allows for separable- variable solutions as 

n(p,(^,C) = <!>(/>, CK"*^. (2.27) 

The evolution of the radial wavefunction $(/), C) for a given azimuthal index m is accordingly 
ruled by the 2D PWE in the radial coordinate p: 



. d d"^ Id m? 
dC, dp^ pdp p2 



$(/5,C) = 0. (2.28) 



In order to apply the above illustrated procedure to this equation, we resort to an appropriate 
representation of the s/(2,]R) ~ sn(l, 1) generators for the symmetry of concern, namely 

As their planar counterpart (|2.4p . the above relate to the operators x^, p^ and x • p, properly 
expressed for a circular cylindrical symmetry. 

The relevant propagator ^pwe(C) •= exp(— i^/C_) amounts to a Hankel-like transform as that 
in rectangular coordinate amounts to the Fresnel transform (j2.3p . One has in fact: 

^{p,0 = [U,^,{Of]{p) = t^ I\i,ip''+P')j^^^2f)f(^p')p'dp', (2.30) 

Jm denoting the Bessel function of the first kind and order m [2^. Of course, #(/3, 0) = f{p). 
The above is a particular form of the (real or complex) radial canonical transform |21j 



$(p,Co) = % Bf{p,a) = ^-^^1 e^(^''"^^'^'h^{^)^p',c^)pdp', (2.31) 

corresponding to the canonical transformation conveyed by the (real or complex) symplectic 
matrix M = (^ ^) . It basically follows from the 2D Collins integral under the circular cylindrical 
symmetry assumption (I2.27p . Hence, it relates the radial wavefunctions of the wavefields (of 
a given azimuthal symmetry) at the input and output planes at Q and Coj between which the 
optical system, described by the ray-matrix M, is conventionally intended to operate. 

As (I2.10p . when A ^ the diffraction integral (I2.3ip can be seen as resulting from a free 
propagation by B/A , followed by a lensing with focal power —C and a scaling by 1/A, so that 

<I.(p,Co) =% BfipXi) = \e'^P'^ (^,C. + I) . (2.32) 

Here ^{v, Qi + B/A) = ^p^^{C,i + B / A)'^ q\{v) represents the radial function composing the prop- 
agated form of the (effective or fictitious) source function u(t;,99, 0) = ^o(^)e*™''^, with which 
the wavefunction u{v,ip,Q) = ^{v,Q)e^'^'^ can be associated and whose azimuthal symmetry, 
accounted for by the factor e**""^ and supposedly preserved by the propagation, is implicitly 
conveyed by the dependence of the free-propagation integral (j2.30p on the azimuthal index m. 
In analogy with (12.17p . the case ^ = yields the transformation 

^(OB)^(p,Ci)= ^ j^ e^'^[^,Q), (2.33) 

involving, as expected, the Hankel transform ^{p/B,Q) of order m of the radial wavefunction 
at the input plane, modulated by the phase factor exp (^2bP^)- 
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As is well known, in fact, the 2D Fourier transform of a function v{p, ip) obeying the factori- 
zation (|2.27|) with respect to the polar coordinates, i.e. v{p,ip) = e*™''' <!'(/?), turns into 

with ^{p) signifying the Hankel transform of ^(p) of order m, which according to the usual 
definition means 

^p) = [Hm^] (p) := / Jm{pp')Hp')p'dp'. (2.34) 

Jo 

As for (|2.2p . several solutions of (j2.28p have been identified, which can be understood as 
arising from eigenstates of definite operators in the algebra generated by {/C+,/C3,/C_} [Tm[T6]. 
Thus, for instance, the eigenfunctions of the operator /C^q = /C3 — (2/^o)^+ evolve into wave- 
functions, which, as the aforementioned Weber-Hermite solutions of (j2.2p . depend on three 
independent parameters and comprise a complex quadratic exponential modulated by the Whit- 
taker first function M^,^ of suitable argument [16]. As remarked, such solutions have been also 
deduced in [21] through an appropriate variable-separation ansatz. 



2.2.1 The Hankel transform and the optical (radial) Appell transformation 

In full analogy with the Fourier transform, the Hankel transform (j2.34p of order m can be given 
an operator representation in terms of the algebra generators /C+ and /C_; namely, 

ii^ = i^+if,-ii^m[!C-+ic+] _ ^2.35) 

As before, the factor i™+^ allows the mathematical transform ()2.34p to be matched to the optical 
transform, conveyed by the Huygens-Hankel integral ()2.3ip for the Fourier matrix F. 

Then, going through the same procedure as before, we can see that the radial Appell trans- 
formation (for a given azimuthal index m) is described by the operators 

A^ ■= e-'^^-Ume'^'^- and A'} := e-'^^-fi^^e'^^- , 

whose dependence on C, will not explicitly displayed. 

They respectively amount to the ABCD matrices (j2.15p and (j2.19p . and hence, by ()2.32p . to 
the (radial) wavefunction transformations (for a given azimuthal index m): 

^(p,C) = ^e^^$(T^,-i). (2.36) 

In accord with (|2.33p . accounting for the case ^ = 0, we recover the Hankel transform-relation 
between the source functions: ^(/3, 0) = $(/9, 0). 

The two possibilities conveyed by (j2.36p actually yield the same expression for the transformed 
wavefunction; this reflects the fact that Am = -^m^ as a consequence of the self-reciprocity of 
the Hankel transform: T-Lm = T~i^. 

In analogy with the "linear" operator A, Am can be regarded as a (back) evolving-location 
Hankel transform operator: Am = T~im{C)- 

Paralleling the cases discussed in Section 12.1.21 we consider, as an Appell pair of solutions 
of (j2.28p . the Bessel beams [l6] and the Bessel-Gauss beams [U]. The former can be interpreted 
as propagated forms of the eigenfunctions Jm{^p), A G M, of the free-Hamiltonian operator /C_: 

Bx,m{pX)=e-'^''^/^Jm{\p), 

the eigenvalue A signifying the transverse component of the wavenumber: A -f-)- fcj^ [TOlfTB]. 
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By the Appell transformation (|2.36p . the -BA,mS turn into the Bessel-Gauss beams: 






5GA,m(p,C) = M^ e*^e'^J„,(A£), 



which in turn arise from the eigenfunctions 5{p — A), A G M, of the dual operator /C_|_ |imil6j. 

As is well known, the Bessel modes are diffractionless (as the Airy beams V'bhj); whilst the 
Bessel-Gauss modes, which as seen evolve from the Hankel transform of the source functions of 
the former, have complementary properties in both the space and spatial frequency domains. 

Evidently, applying the radial fractional Appell transformation (for the given tti), plainly 
understood as 

2" ■= p-KI^-qjOi JC,K.- _ qja / A^ 

one could follow the continuous transformation from the Bessel to the Bessel-Gauss modes at 
any C,. Here, Ti.'^ signifies the fractional Hankel transform, which, resorting to the fractional 
Fourier transform matrix F", by ()2.3ip means 

[n'^^]{p) = [e^(™+i)<^e-*<^[^-+^+l$](p) 

with (j) = aTT/2 G (— vr, vr]. Therefore, as far as cos(j) — (^sin</> ^ 0, A^ yields the transformation 

cos (p — C sm (p \ cos (p — Q sm <p cos (p — (, sm < 

whilst for C = cot (p it acts as 

(1 + C Jsmcp \(l + (, jsmcp / 

2.3 Optical Appell transformation and duality. 

Self-Fourier/self-Hankel and self- Appell wavefunctions 

Since, as reviewed above, the Appell transformation manifests the correspondence between wave- 
functions generated by Fourier or Hankel pairs of functions [9] , it naturally relates to the concept 
of beam duality [351 SS]- However, it does not connect wavefunctions, which are locally dual, 
but connects wavefunctions whose duality traces back to the respective source functions. 

An "Appell transformer" - provided it be implementable - would be so an optical device 
turning a wavefunction into that which one would generate by propagating the Fourier transform 
of the source function of the original wavefunction. In a sense, as the Fourier transformer can 
be understood as a "local dual switch", the "Appell transformer" might be understood as an 
"initial-plane dual switch" . Correspondingly, the fractional Appell transformations individualize 
a family of symmetry transformations parameterized by a continuous parameter, which, filling 
the gap between the "evolution-lines" of a function and its Fourier/Hankel transform, might be 
seen as relating to a sort of "fractional beam duality". 

In the light of the examples considered above, i.e. Airy and Bessel beams, we may visualize 
the problem in terms of two ideal paths, running parallel to each other to depict the (^-lines along 
which the evolution of a given source function and of its Fourier or Hankel transform occur. The 
Appell transformation connects one path to the other at any desired C or better one path to any 



^Note that the comparison should more correctly involve the eigenfunctions of P, i.e. the plane waves e 



iA5 
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other path between the two "extreme" ones, thus allowing us to "have a look" at or to turn the 
wavefunction at hand into a wavefunction whose properties in the space and spatial frequency 
domains are a desired mixture of those of the wavefunctions "lying" on the two extreme paths. 

Evidently, when the source functions are self- Fourier or self-Hankel functions, the two paths 
collapse one into the other. The Appell transformation in fact comes to reproduce at any C 
the direct (or, inverse) self-Fourier/self-Hankel relations obeyed by the source functions. This 
is the case, for instance, of the standard Hermite-Gauss and Laguerre-Gauss modes, which 
respectively arise from the eigenfunctions of the self-dual Lie operators K_ -|-K+ and /C_ -I-/C+; 
as seen above, the latter are the generators of the (both ordinary and fractional) Fourier and 
Hankel transforms. 

In fact, taking into account the explicit expressions for the quoted (normalized) modes, i.e. 

for the standard Hermite-Gauss modes, and 



for the standard Laguerre-Gauss modes, where 

MC) = i + <, 

Hn denotes the Hermite polynomial of degree n and L™ the generalized Laguerre polynomial of 
degree n and order m [29], we can verify that 

l"sFG„(e,c) = i-ir''sHGn{c,c), 

A'^sLGn,m{pX) = {-ir''sLGn,m{pX)- 

The above just reproduce the relations holding between the relevant source functions, which 
are the Hermite-Gauss and Laguerre-Gauss functions sHGn{S,,0) = , e~^ ''^Hn{^) and 

y 2"?i!-v/7r 

sLGn,m{p,0) = \ (n+mY- P^^'^ -^n (p^)' ^^^^ the respective Fourier and Hankel transforms. 

In a sense, the standard Hermite-Gauss and Laguerre-Gauss modes can be considered as self- 
Appell wavefunctions, respectively suitable to a rectangular and circular cylindrical geometry. 

3 Canonical transforms: a short review 

As seen, the interpretation of the Appell transformation within the optical context naturally 
resorts to the Fourier and Hankel transforms [9], which ultimately are integral transforms. 

In general, integral transforms provide a well established and valuable method to solve prob- 
lems in several areas of both physics and applied mathematics. As we know, the roots of the 
method can be traced back to the original work by Oliver Heaviside on the ordinary differential 
equations with constant coefficients occurring in the theory of electric circuits |50(l51j. Due also 
to the availability of large scale computers, the method has then been increasingly extended to 
a wide range of physical and mathematical problems, basically described by partial differential 
equations with assigned boundary and initial conditions [52 1 153 1 [T3]. 

In particular, a special class of integral transforms, named canonical transforms (CT), ap- 
pear widely in optics, in electromagnetism, in classical and quantum mechanics as well as in 
computational and applied mathematics. Their link to the canonical transformations and to 
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the parabolic differential equations as well as to the theory of special functions has been deeply 
analyzed. Fourier, bilateral Laplace, Bargmann, Weierstrass-Gauss transforms as well as Hankel 
and Barut-Girardello transforms are examples of CTs. 

Canonical transformations play a crucial role in classical mechanics. When applying the 
inherent formalism to quantum mechanics, the CTs naturally arise (under specific conditions) 
as associated representations (unitary or not) acting between suitably constructed Hilbert spaces 
of functions on the real or complex domain [SH [55l [56l EZl [58l [59] . 

On the other hand, the initial value problem for evolutionary equations is usually formulated 
in terms of the evolution (or, displacement) operator. The 2D PWE and the (1+1)D SE are 
examples of evolution equations, whose displacement operators realize a unitary mapping of the 
Hilbert space of square integrable functions into itself. The same can be said for the ID HE 
although in this case the relevant evolution operator no longer generates a unitary mapping. It 
is proved that CTs can be realized as evolution operators, generated by second-order differential 
operators through exponentiation to the group by real or complex parameters, and hence as such 
can be associated with evolution equations, ruled by Hamiltonian-like operators (not necessarily 
Hermitian) which are quadratic in the inherent canonically conjugate variables |60| . 

In addition, CTs directly relate to some aspects of the theory of special functions through the 
eigenvalue problem [61]. In fact, special functions, like the Hermite-Gauss, the Laguerre-Gauss 
and the parabolic cylinder functions, are self-reciprocal under some of the aforementioned CTs. 

The Fourier and Hankel transforms are basic examples of CTs. Their role has been enlarged to 
more general contexts by the introduction of the corresponding transforms of fractional order [621 
[63| I64j . which, for instance, in the field of optics [421 H3| [HI I45j . gave rise to a great variety 
of applications, investigations and new formulations in an increasingly enriched optics scenario, 
that in turn stimulated further general analyses of the linear CTs from both theoretical and 
applicative/numerical points of view. Few recent titles are in [65l[66l [671 [68} [69 l [70] . 

The Fourier and Hankel transforms are real (respectively, linear and radial) CTs. In Sec- 
tion 14.21 we will see that the caloric Appell transformation, relevant to the ID HE, naturally 
relates to the bilateral Laplace transform, which in contrast is a special type of complex (linear) 
CTs. 

3.1 Linear canonical transforms 

We recall that real lineaijj CTs realize a unitary mapping of the Hilbert space £'^(R) into itself 
through the integral transform |55[ [56l \E7\ [58] 

/oo 
K^{x,x'Mx')dx'G2.\R), (3.1) 

whose kernel 

K^(x,x') := ^.=^e^(^^"+^^^-2-') (3.2) 

depends on the three linearly independent entries A, B, C of a 2 x 2 unimodular real matrix 
M = (c-^), with the aforestated convention on the square root of a (complex) number and 
the requirement A/B > resulting form the integrability condition. The matrix M belongs 
to the three-parameter real symplectic group Sp{2,M) ~ SL{2,M), and in practice specifies the 
canonical transformation of the involved (classical or quantum) canonically conjugate variableqK. 



The term linear accounts for the integration involving the real line R. 

'''As is well known, classically a canonical transformation is a change of the phase-space variables {q,p) — >■ 
(g'((7,p),p'(q,p)), which preserves the Poisson bracket {q,p} = 1 = {q',p'}. As a straightforward extension of 
the above concept, a quantum canonical transformation is defined as a change of the conjugate non-commuting 
observables g — >■ q'{q,p), p ^ p'{q,p), which preserves the Dirac bracket [q,p\ — i = [q' jp"]. 
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Notably, the inverse transform is given by the same expression with M replaced by the inverse 
M^^: T^-*^ = T , , thus amounting to the kernel K , (x,x') = K* (x',x). 

As said, the Collins integral (J2.1U|) is a real linear CT, the relevant matrix M signifying in 
fact the transformation of the paraxial-ray variables consequent to the propagation through the 
optical system described by M. 

The case B = 0, which optically signifies imaging, yields the geometric transform 

T^^ ^^ : V^ G il\R) ^ ^e^^^V {j) e £' W- (3-3) 

The extension of real CTs to complex CTs is rather involved. For an accurate account, we 
address the reader to the devoted literature [Ml ESI EH EH EEl EH]- Here, we briefly recall that 
complex linear CTs amount to the same integral transform (j3.ip . (j3.2p . but the kernel involves 
the three linearly independent entries of a 2 x 2 unimodular complex matrix M, belonging to 
the six-parameter complex symplectic group Sp{2,C) ~ S'L(2,C); accordingly, the integrability 
condition requires now that '^{A/B) > or B real if A = 0. In fact, when extending the 
concept of canonical transformation from classical to quantum mechanics, it seemed useful to 
extend the transformation from the real to the complex domain. This implies that Hermitian 
operators can be mapped into canonically conjugate, but not necessarily Hermitian, operators. 
Then, contrary to the real CTs, complex CTs no longer represent unitary mappings from £^(]R) 
to £^ (M) ; they transfer 2p (M) into the Bargmann-Hilbert space ^j^ of analytic square-integrable 
functions over the complex plane [5l], completed by a suitably defined scalar product in order 
that the transformed operators have the appropriate hermiticity properties and reproduce in the 
new variable the Schrodinger representation as x and —id/dx. Specifically, the complex CT Tj^ 
associated with the transformation matrix M G Sp{2, C) amounts to the transform pair [57] 



fi^) = ['^mf]i^) = / dx'K {x,x')ip{x'), 

^{x') = I dfiJx)K*Jx,x'Mx), (3.4) 

connecting £^ (M) and 5m ' where 

(^>V'). = / dfi (x)^*{x)i^{x), 
the measure being 

d^l (x) = ./ Ae2^(«^'-2x^*+«*^*') d'Stxd'^x, 

M \ TTV 

witliu = A* D-B*C €C and v = 2'^{B* A) GlR[5llE5lE6lEaEHlE9]. 

A basic example of complex linear CT is provided by the Bargmann transform |54j 

^ ' dx'e~^^'"+^"-^^-^'^'y^ip{x') = {27r)'/^[j^ip]{x), (3.5) 



rV4 

the relevant canonical transformation, signified by the Bargmann matrix 
J_ / 1 -i 



^~V2\-^ 1 

essentially turning the coordinate and momentum operators into the harmonic oscillator raising 
and lowering operators. Accordingly, the Bargmann transform relates the Schrodinger and Fock 
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representations of quantum mechanics; it is in fact applied, for instance, to the coherent-state 
formulation of quantum optics. Evidently, Jg is the space of holomorphic functions on C which 
are square integrable with respect to a Gaussian measure. 

We will specifically deal below with two kinds of complex linear CTs, associated with the 
transformation matrices 

PW=fn "rV ^>0, (3.6) 



1 



and 



G(1M = (,/^ ?)' ^>°- ^^-^^ 

The former produces the Poisson (or, Weierstrass-Gauss) transform 

1 r+oo 

V27rr J-oo 

which may account, for instance, for heat conduction (Section 14. ip . It may be considered as the 
complex counterpart of the Fresnel transform (|2.3|) accounting for paraxial free-propagation, the 
relevant transformation matrix being (j2.1ip . Similarly, it transfers the initial function <p{x) at 
r = into the temperature function 'iIj{x,t) = [Tpip]{x) at subsequent rs. Interestingly, since 
the Poisson transform formally implies the convolution with a pure Gaussian function, it is 
optically implementable by a Fourier transform followed by the propagation through a Gaussian 
aperture followed in turn by an inverse Fourier transform. The parameter 1/r takes the meaning 
as the characteristic width of the involved Gaussian function. 

As to (13. 7p . we see that according to ()3.3p it amounts to a modulation by the Gaussian 
exp{—x'^/2w), which is optically realizable by propagation through a Gaussian aperture, the 
parameter w signifying then the characteristic width of the aperture. Evidently, the real matrix 
L(l//) = [_l/f *! ), amounting to the multiplication by the phase factor exp(— ix^/2/), is the 

real counterpart of ()3.7p : it signifies indeed the propagation through a thin lens. 

As their real counterparts, the above matrices obey the semigroup property: P(ri)P(T2) = 
P(r2)P(Ti) = P(ri+T2) and G(lM)G(l/w;2) = G(l/u;2)G(l/w;i) = G{l/wi + l/w2), and rep- 
resent a dual pair of operators, since P(t) = F~"'^G(r)F and similarly G{l/w) = F~^F[l/w)F. 

An interesting pair of related real/complex CTs is represented by the Fourier and bilateral 
Laplace transform, whose features of interest here will be reviewed below. 

3.1.1 The bilateral Laplace transform 

The bilateral Laplace transform L is well known to transfer an in general complex- valued function 
defined over M into a complex-valued function defined over C by the integral 

/+00 
dx'e-^^V(^'), (3.8) 

-oo 

simply referred to as Laplace transform in the following. 

It can be recovered within the formalism of the complex linear CTs. In fact, according 
to (13. 4p . the unimodular complex matrix 

: ■) . (3.«) 
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yields the integral transform Tl, which is off of L by the factor — i/\/2vr, being 

^ 1 /■+00 1 

[Tl^,] {x) = -j= dx'e— V(x') = -7^ [M (^). (3.10) 



As noted in [57j, since L does not fullfil the aforestated integrabihty condition that B be 
real if A = 0, the relevant transform might be meaningless as the inherent integral might be 
divergent. Indeed, we know that a critical issue in dealing with Laplace transform is convergence, 
since [L(^](x) generally exists only for some values of x, located in the region of convergence 
(determined by ip and 'Six). However, we will not dwell here on such a question, since we 
will not go through the evaluation of the Laplace transform of specific functions, but we will 
exploit its interpretation as a linear CT with the associated transformation matrix (j3.9p and the 
consequent exponential operator representation in terms of the ■su(l, l)-generators K_ and K+. 
In fact, in analogy with the Fourier transform (see (12. 9p ). Tl can be expressed in the form 

fL = e(-/2)[K-K+], (3.11) 

and so interpreted as the non-unitary evolution operator describing the dynamics of a quantum 
repulsive oscillator observed at the purely imaginary "time" r = 111/2 [60^1611 . whereas, as earlier 
noted, the Fourier transform can be understood as the unitary evolution operator describing the 
dynamics of the quantum attractive oscillator observed at the real "time" r = 7r/2. In this 
connection, it may be worth recalling that the operators K_ + K_|_ ex P^ + X^ and K_ — K_|_ ex 
p2 _ y^ span by exponentiation the elliptic and hyperbolic subgroups of the symplectic group. 
Moreover, paralleling the factorizations ()2.13p of J-", we similarly write 

Tl = e^-e-^+e^- = e-^+e^-e-*^+, (3.12) 

since, along with (j2.12p . one also has 

g/3[K_-K+] ^ gtan(/3/2)K_g-sin{/3)K+gtan{/3/2)K_ 



e 



■tan(/3/2)K+gSm(/3)K_g-tan(/3/2)K+^ |^| ^ ^^ 



As is well known, the Fourier transform admits as eigenfunctions the Hermite-Gauss func- 
tions, 



Un{x) = e "'/'//„ (x) 

y 2"n!-^7r 



(that in Section 12.31 have been seen to be the source functions sHG{^, 0) for the standard 
Hermite-Gauss modes), being in fact 

^Un = {-i)"-Un. 

In contrast, the eigenfunctions of Tl involve the Weber-Hermite functions D^/, which for 
integer orders turn into the Uns: Dn{v2x) = ^/ n\y/TTUn{x) p9]. The eigenvalue equation for Tl, 

fL< = Aj/i)ci.;;, 

is in fact solved by 



^t(^) = ^^^|^e-(^-+V2)/4^^^_^^^( i e^3./4^,)^ ^J^) = , 



7r^/2 
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for real values of // [10^ [60l [6T] . The functions $ s provide an eigenfunction basis for £^ (M) as 
well as for the Bargmann-Hilbert space ^^. 

In [nn] a possible definition of fractional Laplace transform has been proposed by applying 
to the fractional Fourier transform matrix F° the same similarity transformation which through 
the ordinary Fourier transform matrix F yields L. In fact, sincqj 

L^g(gi7r/4)FS(e-^'^/^), 

where S is the "dilation" matrix Sie'""/^) = (fl-i^/i) G SL{2,C) - Sp{2,C), one can think of 
reproducing the above scheme with F"^ in place of F. The resulting matrix 



^ ^ ^ ' \t sm (p cos (p J 



has been taken in [66] as the representative matrix of the fractional Laplace transform C"' of 
order a by the definition 

£" = r/2fL„. 
Then, echoing the fractional Fourier transform (see (j2.24p ). C" amounts to the integral transform 



[£«<^](a;) = ./LII^ r"dx'e^(""'=°^'^+"''=°^'^-2""'V(rE'). (3.13) 

V ^v^^ J— oo 

The ordinary transform is recovered with a = 1, which yields C^ = i^'^Ti, = (27ri)~^''^L. 
In analogy with J^", one can prove for >C" the relations 

£"+''=£", i = 0,±l,±2,..., 



the periodicity with respect to the order a allowing us to limit it to the interval (—2, 2], whereas 
the power-like relation implies that C" admits the same eigenfunctions $ as >C (and so Tl): 

It has recently been noted [7l] that the transform £", as proposed in |66| through (j3.13p . 
does not reproduce the fractional Fourier transform (|2.24|) when the real part of the transform 
variable x' = Kx' + z9x' is set to zero. Then, in order to overcome such a limit, a different 
parameter matrix, precisely L'" = Qg^f _i colW ' ^^^ been suggested in the quoted reference 
to define the fractional Laplace transform through the usual scheme (j3.4p pertaining to the 
linear CTs. We have mentioned such a result for completeness' sake, but we will not use the 
proposed matrix L'° in the forthcoming analysis. Actually, also L" will marginally be used 
below. 



*It may be worth noting that it is also L = S(e"^''^)F = FS(e~"^/^), and hence Tl relates to Tf with the 
relevant integral involving the imaginary axis. As remarked in [57], in fact, multiplication of the transformation 
matrix M on the left (on the right) by the "dilation" matrix S(e'") (S(e~'")) is a tool to obtain transforms 
involving line integrals along a path tilted by the phase a. 
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3.1.2 Some about the Laplace transform representative matrix 

The representative matrix of the Laplace transform, be it L or L", if we are dealing with the 
fractional transform (I3.13p . belongs to the set of complex unimodular matrices as 

M=(_-^ o). P-") 

with A, B, C, D real, and AD — B C = 1, which indeed form a subgroup since the product of 
two matrices of this type yields a matrix of the same form. 

Complex canonical transformations, conveyed by matrices like ()3.14p . and their nonunitary 
representations have a number of interesting applications in physics. They arise, for instance, 
in the clustering theory of nuclei |59j as well as in the study of the accidental degeneracy [72] in 
both a 2D anisotropic oscillator, whose frequencies in the two directions have a rational ratio, 
and a 2D isotropic oscillator, constrained to move in a sector of angle Tr/q, q integer. 

Also, the matrices (|3.6p and (|3.7p . representing optical processes like the convolution by 
a Gaussian function and the propagation through a Gaussian aperture, are of the type (13.14p . 

Specifically, the matrices (j3.14p produce the integral transform 

[f ^] (x) = ^4= / e^(^^''^^^'-'^^yix')dx'. (3.15) 

As the Huygens-Fresnel transform (j2.10p . it can be recast in a form similar to (|2.16p . In fact, 
(j2.16p reflects the factorization a la Wei-Norman of real symplectic matrices |28j , 

'A B\ _ ( 1 0\ (A \ (I B/A\ 
^C d) - [C/A l) [o l/Aj [o 1 ; ' 

in terms of a dilation matrix sandwiched between two real triangular matrices, respectively 
lower-left and upper-right, which optically account for lensing and free-propagation. 
A similar factorization of matrices of the type (|3.14p is allowed as 

A iB\ ( \ 0\ M \ (\ -ir' 



-ic D) \i/w ly vo iMy vo 1 y ' ^^'^^^ 

with the real parameters r and w given by 
_ _B_ 1__^ 

Thus, by (j3.16p the transform (j3.15p can alternatively be understood as 

where 99(x,ro — B/A) is to be intended as the "Poisson" -transformed (i.e. in a sense, "evolved" 
with respect to the effective or fictitious parameter tq) form of the function (/j(x), namely 

1 r+oo ^ 

^(x,ro + t) = — == / e-2^^''-'^'^\o(x)dx = ^^ipo]{x). 
V^vrr J-oo 

with >fQ{x) = ip{x, To). 

As for real matrices, the factorization (|3.16p holds for A ^ 0. If ^4 = 0, the general expres- 
sion ()3.15p yields 

which involves the Laplace transform (j3.8p of the function. 

Of course, the interpretation (|3.17p of (|3.15p can as well be obtained by directly manipulating 
the expression of the transform. 
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3.2 Radial canonical transforms 

Besides the linear transforms, radial transforms involving (square-integrable) functions defined 
on the positive half-line M~^ can as well be associated with canonical transformations. 

As theorized in [58], the radial CT T^^'"*^ generated by the canonical transformation of the n- 
vectors x = {xj} and p = {pj}, j = 1, . . . , n, which, preserving the Dirac brackets [xj,pk] = iSj^k-, 
is identified by a 2n x 2n unimodular complex matrix M = (cd)' whose entries are (complex) 
multiplies of the n x n unit matrix such that AD — BC = 1, amounts to the transform pair 

^(0 = [Ti:j'™V](0 = / dr'K(;^'™)(r,r')r'"-V(r'), 

Jo 

^{r')=f d^,^^^\r)K'^^'^\ryr^{r), (3.19) 

the kernel being 

By (j3.19p . functions ^{r') in the Hilbert space £^(M"'") with scalar product 

(¥','0)£2(K+) = / r"'-^ip*{r')%l){r')dr', 
Jo 

are transformed into functions ^(r), belonging to the space d^'"^' of analytic functions of the 

complex variable r restricted to the region C"*", for which ■^ip{r) = 0. Moreover, the space 

'^M ^^ equipped with the scalar product 

('?,V^)(_, = / rfA^^''")(r)^*(r)V^(r), 
^m Jc+ "^ 

the measure being 

d^{n,m)r^ ^ A g 2^ («^^' +«*'•*') (rr*)"/2i^„ /2+„,i ( —] d^r d^r, 



where K^ denotes the MacDonald function of order v [29], and, as before, u = A*D — B*C ^ C 
and V = 2Q{B*A) G M \57\ I58j . It ensures that the transformed operators have the appropriate 
hermiticity properties and reproduce in the new variable the Schrodinger representation of the 
operators x^ = r^, x • p = —ir-^ and p^ = —{-ifi + ^^^ + ^)) with the real parameter 
A = —m(rn + n — 2), m = 0, 1, 2, . . . , addressing the eigenvalue of the angular momentum L^ = 
2 Yl LjkLjk, Ljk = XjPk — XkPj, which is invariant under the unimodular transformation M |58j . 

In particular, with n = 2 and M resorting to the entries of the Fourier transform matrix F, 
one basically obtains the Hankel transform (j2.34p . whereas the fractional Fourier transform 
matrix F"^ yields through (|3.2U|) the fractional Hankel transform (|2.37p . 

The Barut-Girardello transform G„,m represents an interesting example of complex ra- 
dial CT. Resorting to the Bargmann matrix B, it signifies |73j 

r+oo 

[Gn,mV] {r) = V2 rfr'(rr')i-"/2e-i(^'+^")4/2+^_i(v^rr')</.(r') 

Jo 

= i'^/2+'»-i [f |^"''")c^] (r), (3.21) 

Ii, denoting the modified Bessel function of the first kind of order u: Jy{ix) = i^Iy{x) |29j . 
It was introduced in developing the formalism for coherent states associated with Lie algebras 
of non-compact groups, like, in particular, the semisimple Lie algebra so(2, 1) ~ sZ(2,M). As 
a direct generalization of the notion of coherent states associated with the Heisenberg algebra, 
such generalized coherent states were introduced as eigenstates of the lowering operator of the 
aforesaid algebra in the relative discrete representations D^{j), j = —1/2, —1, —3/2, . . . [73] . 
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3.3 Essentials of Hankel-type transforms 

The first and second Hankel-type transforms of Bessel order v > —1/2, depending on an arbi- 
trary real parameter u' , are respectively defined by [74l [75] 



[Kyf](y^^k^.Ay)--=y'''"'' / {xy)-''M^y)f{x)dx, 



+ 00 



[^,..j](y) = /,..,(y):= / x'+''''ixyy^'Mxy)f{x)dx, (3.22) 

Jo 

the function f{x) being assumed to belong to the space £^(M"'"). 
It can be seen that they relate to each other according to 

7x , ^ 7x„ , , rL , ^ ri, , . 

1,1^, ly' 2,1/, I/' 2,1/, I/' l,!/,!/' 

Also, as a consequence of the well-known orthogonality relation of the Bessel functions, both 
transforms are self-reciprocal: 

l,v,v' 1,1/, I/" 2,!/,!/' 2, 1/, I/" 

and separately obey the Parseval relations |76j 

hoo r+oo ^ 

x~^~^^' nx)9{x)dx = x-^-^'^'r ,{x)g^^Jx)dx, 

Jo '"'" 

x'+'^'r{x)g{x)dx= / xi+2A«/;_,(x)5,_,(rr)dx, 
along with the mixed relation 

-OO r+oo 



f*{x)g{x)dx= / f*^^Xx)g^^^,{x)dx. 

U "^ u 

In [77] the corresponding fractional order transforms have been introduced. 
Notably, one can prove the identities (711 EH] 

[^.,..'^!„v/](y) = -y^[K.,J]iy)^ [KyK.'f](y^ = V[^,,.,.,/](y), (3.23) 

involving the Bessel-type differential operator B , and its adjoint B' , : 
"'"' dx^ X dx x"^ dx dx 



St , := -^ - (1 + 2u')-^ + [(^' + 1)2 -v']\- .^'- A.2.+1 A.-'— 1 



dx^ X dx ^ -' x^ 9x 9x 

Relations (|3.23p can be used, for instance, to write the solution of the evolutionary equation 
[711[75] 

K^h{x,T) = Bl^Xx,T), (3.24) 

satisfying the initial condition /i(x, 0) = /(x), in the transform conjugate y-space as h , (y, r) = 
^-{T/K.)y f^^^^y^ for arbitrary values of k. Then, transforming back to the x-space, we obtain 

r+oo 

h{x,T) = -x'^^'"' j {xy)-''e-^^^"+y"h,{^xy)f{y)dy, 
under the condition that |arg(r/K)| < 7r/4, which for both r and k real turns into t/k > 0. 
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As noted in [77], on account of the symbolic solution h{x,T) = exp[{T/K)B^ ,]f{x) of p.24p . 
the above yields an explicit functional representation of the exponential operator expf/3-B^ 1: 

[/""lyfjix) = ±x'+^''' l^^{xy)-^'e-M^'^y"h,{^)f{y)dy. (3.25) 

For /3 = i/2 it eventually conveys a representation of the transform T-i_^ ^ ^, as a symmetric product 
of exponentials of the su(l, 1) algebra generators, similar to (j2.35p for the Hankel transforms: 

This is not surprising since the operators 

K^):=l.^ K^^=4^l. K^^^=-§(x|:-.'), (3.27) 

constitute a non self-adjoint one- variable realization of the su{l, 1) algebra generators obeying, 
as it should be, the inherent commutation relations shown in ()2.4p . 
The adjoint operators 

pertain to the second Hankel-type transform T-L,^ ^ ^, , which is in fact amenable for a factored 
form representation analogous to (I3.26P in terms of them [77j . 

With v' = —1 and v = m the operators ()3.27p turn into ()2.29p . pertaining to the Hankel 
transform T-Lrn of order m, considered in Section [2.21 in connection with the propagation problem 
in circular cylindrical coordinates; indeed, T-L^ ^^/ ~^ ^m for the identified values of u and u' . 

In general, as noted in [77], for suitable values of u' and u the transforms T-i^ , and 11^ , 
can be framed within the formalism of the radial CTs, as theorized in [58] and shortly reviewed 
in Section [3. 2[ In fact, if the 2n x 2n transformation matrix resorts to the entries of the Fourier 
matrix F, the kernel p.20p becomes 

K^n,H(,,/) = (-.)'"+"/2(rrOi-"/V„/2+^_i(rrO, 

the resulting transform being seen as the radial part of a n-dimensional linear CT, specifically 
representing a 7r/2-rotation for each pair of the canonically conjugate operators in the relevant 
n-component position and momentum operator vectors. The kernel K^"'™' resembles that of 
the transforms ()3.22p for suitable settings of the parameters v and v' in terms of the inherent 
dimension n and eigenvalue A. In fact, one respectively has n = —2v' and n = 2(1 + z/'), which 
so are meaningful for v' integer or half-integer (negative or positive), whereas correspondingly 
771 = 1/ + I/' + 1 and m = V — v' , thus yielding \= {v' + 1)^ — z/^ and X = v'^ — v^ . 

Disregarding the specific link of the canonical transforms to the canonical transformations as 
briefly reviewed above, we rewrite the integral transforms ()3.22p through the formal relations 



W.Mir) = ^'-^^ / K(-2^''^+^'+i)(r,r')r-i-2^V(/)dr', 
Jo 

K(2(i+.'),-^')(,,/)/i+2.'^(^/)^^/^ (3.29) 



which, contrary to the general account (I3.19p . should be intended as transforming functions of 
a real variable into functions of a real variable. 
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In virtue of the analysis presented in Section [2l one may guess that the above might be of 
relevance in connection with the Appell transformation for a "radial PWE" in a higher-dimension 
domain. 

In Section 14.31 we will prove that the "Laplace-like" versions of (|3.29|) relate to the Appell 
transformation for the radial HE, which is seen to connect temperature functions evolving from 
functions linked through the aforesaid transforms. 

4 Caloric Appell transformation: 

Laplace and Hankel-type transforms 

In virtue of the well-known analogy between the ID HE and the 2D linear PWE, we will properly 
reproduce within the context of the ID HE the interpretation of the Appell transformation, 
as elaborated in Section 12.1.21 in connection with the optical propagation. We will find that 
the well-known caloric Appell transformation connects temperature functions resulting from 
initial conditions linked by a Laplace-like similarity transformation. Thus, as the optical Appell 
transformation is understood to manifest the action of the evolving-location Fourier transform 
operator on wavefunctions, likewise the caloric Appell transformation will be seen to result 
from the action of the evolving-time Laplace-like transform operator on temperature functions. 
A similar relation will then be shown to hold also for the Appell transformation for the radial 
HE, accordingly resorting to suitably defined Hankel-type transforms, explicable as "radial- 
Laplace" -type transforms. 

4.1 ID heat equation 

The ID HEJHl 

2^«(^'i) = ^"(^'i) (4-1) 

has been the object of extensive studies; few titles from the wide bibliography are in [H [21 [3l 
HI El m [7]. See also [HI [131 [H]. The temperature or caloric function u{x,t) is usually assumed 
to be of class C^ with respect to t and of class C^ with respect to x across the region V in 
the x,t-plane where (j4.ip holds [2j. In view of the analysis we will perform, u{x, t) will be here 
supposed to be of class C°° with respect to the variables x, t, which will be considered as real 
unitless variables. 

As is well known, solutions to ()4.ip can be obtained as the Poisson transform of the initial 
conditions, that is [2] 

1 r'^^ I _ ■)2 
/(x,t) = [U„,{t)fo]{x) := -= / e-^^^fo{y)dy, (4.2) 



with /o(x) = /(x,0), under the general statement that fo{x) be Lebesgue integrable in every 
finite interval. The integral kernel reproduces the fundamental solution of (|4.ip : 

1 „2 



S{x,t) := ^^e"2F, t > 0. 



* Actually, the standard form of the ID HE is [2] 



dt ^ ' ' dx 



7J7'^(^'*) = 7n5"(^''*)- 



Here, in order to favor the comparison with the 2D PWE, we will work with (|4.ip . which amounts to the trivial 
variable scaling: t — >■ f/2 between the "time" in the standard form and that in (|4.1|l . 
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Resorting to the review in Section 13. H we see that U^^ = Tp , with the replacement r — ?• t in 
the relevant canonical transformation matrix P. Evidently, the Poisson transform (|4.2|) corre- 
sponds to the Fresnel transform ()2.3p by the variable interchange i o i^, which is the basis of 
the analogy between the ID HE and the 2D PWE. It conveys an integral transform represen- 
tation for the "evolution" operator ZY^^ (t) associated with (j4.ip , which can as well be given the 
exponential operator representation 



i d 



U^,{t) := e^l^ = e-'^- , (4.3) 

in analogy with that of the optical propagator ZYp^g(C). As in that case, we consider the two 
representations as fully equivalent, disregarding indeed the question of specifying the domain of 
functions over which such an equivalence effectively holds [2]. 

We have expressed U^^{t) in terms of the operator K_, thus signalizing that in connection 
with the ID HE we will resort to the same basis operators (j2.4p as for the 2D PWE, although 
linear combinations of them by in general complex coefficients will more properly be consid- 
ered. This will amount to deal with complex canonical transforms. In fact, on discussing the 
caloric Appell transformation we will be led to consider the Laplace transform represented by 
the complex symplectic matrix L (Section 13.1. 1|) . whereas the afore-discussed optical Appell 
transformation has led to the Fourier transform, which in turn resorts to the real symplectic 
matrix F (Section [2X2]). 

4.2 Some about the heat polynomials and the associated functions 

The heat polynomials Vn{x,t) are defined as the polynomial solutions of (|4.ip . which realize the 
power series expansion of the simple exponential solution [2], i.e. 

e^-+^"^ = y^^vn{x,t), (4.4) 

n=0 

with X arbitrary parameter. They explicitly write as 



[n/2] ■ 

Vn{x,t) := n\ y —^, TT-^"^^-"' 

^„2^j!(n-2j)! 



X 



which reveals their connection with the Hermite polynomials as Vn{x^ t) = {—t/2)^''^Hn{x/\/—2t). 
Interestingly, the VnS solve ()4.ip subject to monomial initial distributions Vn{x,0) = x"', 
being [2] 

1 roo ,_ ,2 

Vn{x,t) = / e '^^y'^dy, < t < oo. (4.5) 



In addition, they obey the second-order differential equations 
(f d \ 

which determine the x-dependence of the v^s, whose complete expressions follow then from the 
inherent initial conditions u„(x,0) = x". 

The associated functions Wn{x,t), which are solutions of the HE as well, are obtained from 
the VnS by the caloric Appell transformation [Il[2]: 

Wn{x,t) ■.= S{x,t)Vn{j,-j\ = f" S{x,t)Vn{x, -t). (4.7) 
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They admit the generating function relation 

1 (^-x)^ ■^ x^' 



Z / an} 



n=0 

which parahels (|4.4p . as weh as, in analogy with (|4.5p . the integral representation [2] 
1 /"^^ 

^71" J-oo 

Finally, the WnS obey the differential equation 

t-—,+x— + n + l\wn{x,t) = Q. (4.8) 

We address the reader to the bibliography [2l O IH [5] for a more complete account of the 
topic. 

In [9] the optical analogs of the heat polynomials and the associated functions have been 
deduced, which have suggested the interpretation of the optical Appell transformation in terms 
of Fourier-related source functions, as illustrated in Section [2.1.21 

4.3 Laplace transform and Appell transformation for the ID heat equation 

As earlier remarked, the suitable Lie-algebra context should be now that of the six-parameter 
algebra sZ(2,C); accordingly, taking account also of the Weyl algebra, we should deal with the 
semigroup WSL{2, C). However, in our analysis we will make use of CTs of the type (|3.15p (as- 
sociated with matrices of the form (I3.14|) ). Hence, we can limit ourselves to deal with SL{2, C), 
even better with the subgroup C SL{2, C) formed by the matrices ()3.14p . that amounts to a sub- 
algebra C s/(2,C), spanned by the basis operators |K+,K3,K_} through linear combinations 
with coefficients respectively of the type {za,7, i/3} with a, /3, 7 real. 

We consider for every operator K in the aforesaid subalgebra C sl{2,C) the relevant "dual" 

operator K = TlKT£ , obtained from the former by a similarity transformation involving the 

Laplace-like operator Tl. Due to (jS.lip or ()3.12p . K belongs to the same subalgebra as K. 

By retracing the eigenstate-based procedure [10l[T2] applied to the 2D PWE in Section [2.1.2l 
we consider the eigenvalue problems for K, as represented by the equivalent eigenvalue equations 

K/a(x) = A/a(x), (4.9) 

or 

K{t)vx{x,t) = Xvx{x,t), (4.10) 

and for K, as similarly described by 

Kgxix) = Xgx{x), (4.11) 

or 

K{t)wx{x,t) = \wx{x,t). (4.12) 

Here, the evolving-time operators explicitly mean 

K(f) ■.= u,,{t)ku;;^\t) = e-'^-Ke'^-. 
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Just as for the PWE, the above equations signify that we are looking for temperature functions 
v{x,t) and w{x,t) corresponding to initial conditions linked by the integral transform (|3.10p . 
i.e. v{x,Q) = f{x) and w{x,Qi) = g{x) = [Tl/](2;), with in particular f{x) being identified as an 
eigenstate of a given operator in the subalgebra of concern. Hence, as v\{x,t) = e~*^- fx{x), 
the w\s are in turn obtained from the v\s by the "instantaneous" transformation 

wx{x,t) = e-'^-gx{x) = e-'^-JLe'^-vx{x,t) = Avx{x,t). 

The inherent symmetry operator A, identified by the same symbol as the optical Appell operator, 

l:=e-*^-TLe*^-=TL(t), 

signifies now the evolving-time Laplace operator. Its explicit dependence on t will be omitted. 
Contrary to the optical operator, the caloric Appell operator individualizes the complex 
canonical transformation matrix 

fl -it\ fO A A it\ ft iil + t"^^ 

Being it of the type (j3.14p . we can apply the recipe (|3.17p with tq = t, thus inferring the associa- 
ted functions wxs from the vxs in accord to the Appell rule (14. 7p (apart from the factor l/v2vr): 

wx{x,t) := Avx{x,t) = —=e^^vx [t^-j 

By ()3.18p one recovers the primary relation between the "source" functions, i.e. g{x) = [Ti,f]{x). 

Notably, disregarding the formal complexity of the inverse transform and naively using the 
inverse matrix M2;_i = (~_. ^*^ ^ '^, one obtains the correct rule to recover the vs from the ws. 

As a basic example, we may consider the heat polynomials and the associated functions. 
As signalized by ()4.5p , the former are caloric functions arising from the diffusion of the mono- 
mials x"s. Therefore, the eigenvalue equations (|4.9p and ()4.10p must be specialized for the scale 
transform generator K3: 

K3X" = Ax", 

with A = —i{n + 1/2), n = 0, 1, 2, . . . , and for the relevant evolving-time operator K3(t): 

K3{t)vn{x,t) = [K3 + itK^]vn{x,t) = Av„(x,i), 

eventually resulting into (14. 6p . obeyed by the VnS. 

As to the dual operator K3 = TlKsT^ = — K3, the eigenvalue equation (|4.1ip writes as 

K35A = -Afi-A, 
whereas (j4.12p for the evolving-time operator K3(t) = — K3 — itK_ specializes as 

[K3 + itk_]wn{x,t) = 2 ('^+ 2 ) ^'"(^'*^' 

which, as expected, yields ()4.8p obeyed by the associated functions Wn- 

Thus, the caloric Appell transformation of the heat polynomials (which, as seen, are asso- 
ciated with the eigenvalue problem for K3 with A = —i{n + 1/2), n = 0, 1, 2, ... ) yields the 
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temperature functions associated with the eigenvalue problem for the Laplace dual operator 
K3 = TlK3T£ , or equivalently with the eigenvalue problem for K3 with opposite eigenvalues. 

Notably, since the Fourier and Laplace dual operators of K3 are the same: J-K3T~^ = 
TlK3T£ = — K3, as discussed in [9], the Appell transformation of the heat polynomials can as 
well be associated with a Fourier-similarity relation, whilst in general, as seen above, one must 
resort to a Laplace-similarity relation when dealing with generic temperature functions. 

In order to complete the correspondence with the 2D PWE, we can say that a fractional 
caloric Appell transformation may also be introduced as the evolving-time fractional Laplace 
transform operator 

i" := e-*^-fL«e*^- =fL-(i), 
with Tlq given as discussed in Section [3.1.11 Evidently, it identifies the transformation matrix 
'1 —it\ fcos(j) ism(j)\ /l it\ _ f cos (p + t sin (p i(l -|- t^) sin( 



-^" \0 1 J \i sin (f) cos (p J \0 1 J \ i sin (f) cos (p — t sin < 

and hence maps temperature functions into temperatures functions according to 

I ^\ 7a , ^^ 1 ~rr, — TJ ■ .. .r^ / X t cost/)- sin ( 



^cos (p + t sin (f) \ cos 4> + t sin <j) ' cos (j) + t sin 

4.4 Hankel-type transform and Appell transformation 
for the radial heat equation 

The radial HE is signified by 



d_ 

'dt' 



2—u{r,t) = D^u{r,t), (4.13) 



where /i is an arbitrary real parameter and D^ denotes the differential operator D^ := #^ -|- 
^^^ [Z3- When /i = n, a positive integer, D^ becomes the Laplacian operator in the radial 
coordinate appropriate to a n-dimensional Euclidean space: Dn = -^ -|- • • • -|- #t- As before, 

we consider both r G [0, -|-oo) and t as real unitless variables; also, as in (j4.ip . for computational 
convenience we have added the factor 2 with respect to the equation appearing in the literature, 
thus implying the "time" scaling: t — >■ t/2 with respect to that equation. 

For /i > 1, ()4.13p admits as solutions the radial heat polynomials Rn,^{r,t), n = 0, 1,2. . . , 
and their Appell transforms Rn,^{r,t), explicitly given by [78] 

i2„,^(r,t):=2"n!i"L(^/2-i('-^ 



r 1 

where Sf^{r,t) denotes the fundamental solution of (j4.13p : 

1 -2 

(27rt)'^ 



Rn^r, t) := S^{r, t)i?„,^ ( 7, -7 ) = t'^^'S^ir, t)i?„,^(r, -t), (4.14) 



Sf^i^^i) ■■= 7:r-7JU2^ "• 



With /i = 2 we can recognize a formal correspondence between (j4.13p and the 2D PWE (j2.28p 
with m = 0, which so rules the evolution of the radial wavefunction for a circularly symmetric 
wavefield. Indeed, just as solutions to the 2D radial PWE follow from the Hankel-Poisson 
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integral (j2.30p under definite initial conditions, likewise solutions to ()4.13p may be obtained 
from assigned initial conditions u{r, 0) = (p{r) by the integral transform 

r+oo 
u{r,t) = Kf,{r,r';t)ip{r')dr', (4.15) 

^0 



with 



K,(r,r';t) := iri-^/VW2e-^(^-^+^'^)j^/,_^ ^!:^^ 



which yields, for instance, the heat polynomials Rn,ii{f,t) with (/j(r) = r^"'. 

More precisely, we may recognize a correspondence with the Hankel-type transforms, dis- 
cussed in Section 13.31 which, as there noted, for specific values of the relevant parameters yield 
the conventional Hankel transform, or more general (real and complex) radial transforms. In 
fact, on account of the above clarified link between the ID HE and the Laplace transform, it 
is natural to consider the kernel one would obtain from (j3.20p with the entries of the Laplace 
matrix L, 

K[".™)(r,r') = (-)"^+"/2(rr')i-"/2/„/2+^-i(rr'), 

which then, through (j3.19p . should yield a sort of "radial" Laplace transform. More in general, 
using it in (j3.19p with the aforestated correspondences {n,m) f-)- (i/, i/'), one would obtain in 
a sense the complex versions of the transforms 7i , and % , , i.e. 

Jo 

f+oo 

[^...Mi^) ■■= i-)'^" / r''^"''{rr')-'''lu{rr'Mr')dr'. (4.16) 

Jo 

Clearly, they relate to the Laplace transform just as the Barut-Girardello transform (I3.2ip 
relates to the Bargmann transform ()3.5p . In this connection, we mention that in [77], reproducing 
the structure of the transforms ()3.22p . two analogous Barut-Girardello- type transforms have 
been introduced along with the relevant fractional-order versions. 

It is easy to see that relation (j3.25p yields for C^ ^ ^, the operator representation 

in terms of the operators (j3.27p : of course, a similar representation can be obtained for the 
adjoint transform L^ ^ ^, in terms of the operators (I3.28p . 

The correspondence between (j4.17p and (j3.1ip is evident. 

It is also evident that relation (|3.25p (as that involving the adjoint operator B ,) can be 
made to reproduce the general expression ()4.15p for the solutions of the radial HE by setting 
V = ii/2 — 1 and v' = — /i/2 (as correspondingly, u = u' = /i/2 — 1). In fact, writing, for instance: 

^RHE(i):=e-*^-\ (4.18) 

with K = — ii?^, , = —^Dn, we see that, by 

u{r,t)= [l(^^,it)^]{r), 
one recovers (j4.15p on account of (j3.25p . 
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Then, just reproducing the considerations previously developed in connection with the 2D 
PWE and the ID HE, in the light of the correspondence of (jiTTl) to (fXTT]) and of (jIT^ to (fOj) . 
we can argue that the Appell transformation of solutions of the radial HE arises from the action 
of the evolving-time "radial-Laplace" -type transform operator L^ n-i- ni^) °^ temperature 
functions. Explicitly, the relevant Appell operator ^^ is to be intended as 

A — e-*^-'£ e*^-' - C it) f4 191 



iK'^) 



which amounts to the parallel "heat conduction" problems: u{r,t) = e - ^{f) and w{r,t) 



e-*K- V(r), with V(r) = Aw2-i,-./2'^(^)- 

The direct evaluation or the use of the matrix-based procedure leads to the well-known recipe 
(see gH])): 

A^u(r, t) = t^^'/^e^^u ( -, -- j oc 5^(r, t)u ( -, - j 

thus confirming the interpretation (I4.19p . 

Needless to say, as before a fractional Appell transformation can as well be introduced, by 
exploiting the fractional versions of the transforms ()4.16p , inspired by those of the Hankel-type 
transforms |77| . 

5 Concluding notes 

After reviewing the interpretation of the optical Appell transformation [9] as conveying the corre- 
spondence between solutions of the 2D PWE, which evolve from Fourier or Hankel-related source 
functions (according to whether the propagation problem concerns a rectangular or a circular 
cylindrical geometry), a similar interpretation has been shown to hold for the caloric Appell 
transformation. It manifests indeed the correspondence between temperature functions arising 
from "source functions", which are related by a Laplace or a "radial-Laplace" -type transform 
according to whether the ID HE or the radial HE is concerned. The optical Appell transforma- 
tion is a symmetry transformation for the 2D PWE, which can be understood as an evolving- 
location Fourier or Hankel transform operator. Similarly, the caloric Appell transformation is 
a symmetry transformation for the HE, which can be regarded as an evolving-time Laplace or 
"radial-Laplace" -type transform operator. 

Also, resorting to the fractional versions of both the Fourier/Hankel and Laplace/Hankel-type 
transforms, one can introduce a fractional Appell transformation in relation to the PWE [9] as 
well as to the HE, thus displaying a family of symmetry transformations for both equations, 
parameterized by a continuous parameter. 

The analysis has resorted to the Lie-algebra based method, as originally developed in a series 
of seminal papers by Kalnins, Miller and Boyer |10t \TT\ [T2] in connection with the time-depen- 
dent SE. Also, it echoes the analysis mastered in jl7| . aimed at characterizing transformations 
between wavefunctions in terms of transformations of the respective source functions. 

As earlier noted, the Appell transformation has relevance in other contexts as well. It would 
be interesting to establish whether it is amenable for an analogous interpretation even when 
framed in other contexts, like, for instance, that of the Kolmogorov equation, considered in [8]. 

We conclude saying that other results concerned with the HE can be revisited within the opti- 
cal context. For instance, the property proven in [6], according to which the Appell transforma- 
tion is essentially the only transformation mapping solutions of the (in general, n-dimensional) 
HE into solutions, applies to the optical context as well. Thus, every transformation mapping 
wavefunctions into wavefunctions can be understood as composed of Appell transformations. 
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scalings and shifts of all the variables. Leutwiler's result has been restated in |18j by resorting 
to the ray-matrix formalism of paraxial optics, and accordingly exploiting suitable matrix fac- 
torization tools, like, for instance, the possibility of realizing any 2x2 optical ray-matrix, which 
ultimately identifies a symmetry operator for the 2D PWE, by an appropriate sequence of free- 
sections and lenses, that then can be rearranged to yield the Appell transformation matrix (j2.15p . 

Since scalings and shifts of the inherent variables do not change the functional form of the 
functions acted on by them, we can say that the world of the solutions of the 2D PWE can be 
understood as formed by two classes of functions, which are just mapped one onto the other by 
the Appell transform. Further solutions can be reached by scalings and/or shifts of both the 
variables ^ and C; without changing evidently the functional form of the original wavefunction. 

This is not surprising, since, as noted in [18], the symmetry algebra pertaining to the 
PWE (j2.2p arises from Fourier-similarity related operators, being in fact J^XJ^~^ = — P, 
^P^-i = X, and accordingly TK^T'^ = K^, TK^T'^ = -K3. 

Ultimately, the Appell transformation is a manifestation of such a property. 
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